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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 

RESEARCH  MEMORANDUM 

ON  THE  RANGE  OF  APPLICABILITY  OF 

THE  TRANSONIC  AREA  RULE 

By  John  R.  Spreiter 

SUMMARY 


Some  insight  into  the  range  of  applicability  of  the  transonic  area 
rule  has  been  gained  by  comparison  with  the  appropriate  similarity  rule 
of  transonic  flow  theory  and  with  available  experimental  data  for  a  large 
family  of  rectangular  wings  having  NACA  63AXXX  profiles.   In  spite  of  the 
small,  number  of  geometric  variables  available  for  such  a  family,  the 
range  is  sufficient  that  cases  both  compatible  and  incompatible  with  the 
area  rule  are  included. 


INTRODUCTION 


A  great  deal  of  effort  is  presently  being  expended  in  correlating 
the  zero-lift  drag  rise  of  wing-body  combinations  on  the  basis  of  their 
streamwise  distribution  of  cross-section  area.   This  work  is  based  on  the 
discovery  and  generalization  announced  by  Whitcomb  in  reference  1  that 
"near  the  speed  of  sound,  the  zero-lift  drag  rise  of  thin  low-aspect-ratio 
wing-body  combinations  is  primarily  dependent  on  the  axial  distribution 
of  cross-sectional  area  normal  to  the  air  stream."  It  is  further  conjec- 
tured in  reference  1  that  this  concept,  known  as  the  transonic  area  rule, 
is  valid  for  wings  with  moderate  twist  and  camber.   Since  an  accurate  pre- 
diction of  drag  is  of  vital  importance  to  the  designer,  and  since  the  use 
of  such  a  simple  rule  is  appealing,  it  is  a  matter  of  great  and  immediate 
concern  to  investigate  the  applicability  of  the  transonic  area  rule  to  the 
widest  possible  variety  of  shapes  of  aerodynamic  interest. 

The  experimental  data  contained  in  reference  1  and  many  subsequent 
papers  have  shown  that  this  simple  rule  is  often  remarkably  successful  for 
a  wide  variety  of  shapes  ranging  in  complexity  from  simple  bodies  of  revo- 
lution to  models  of  complete  airplanes.   Furthermore,  important  reductions 
in  the  transonic  drag  of  wing-body  combinations  have  been  realized  by 
indenting  the  body  so  that  the  axial  distribution  of  cross-section  area 
corresponds  to  that  of  smooth  bodies  of  revolution  having  low  drag  at 
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supersonic  speeds.   On  the  other  hand,  it  is  important  not  to  overlook 
the  fact  that  there  are  a  number  of  test  results  on  equivalent  bodies  for 
which  the  correlation  of  drag  rise  by  the  area  rule  is  unsatisfactory. 
Inasmuch  as  the  models  tested  are  generally  of  complex  geometry,  and  only 
the  original  model  and  an  equivalent  body  of  revolution  are  tested,  it  is 
difficult  to  ascertain  whether  these  discrepancies  are  attributable  to 
viscous  phenomena  or  to  the  fact  that  the  drag  rise  may  depend  on  other 
geometric  parameters  than  the  axial  distribution  of  cross-section  area. 

It  is  the  purpose  of  this  note  to  examine  in  further  detail  the 
applicability  of  the  area  rule.   Despite  the  fact  that  most  of  the  emphasis 
in  the  tests  relative  to  the  area  rule  has  been  on  wing -body  combinations, 
there  exists  such  a  scarcity  of  experimental  data  of  a  sufficiently  sys- 
tematic type  that  the  present  discussion  will  be  confined  to  rectangular 
wings  without  bodies.   Transonic  drag  data  are  available  from  bump  tests 
in  the  Ames  l6-foot  high-speed  wind  tunnel  for  a  large  family  of  rectangu- 
lar wings  having  NACA  63AXXX  sections,  aspect  ratios  varying  from  -0.5  to 
6.0,  thickness  ratios  from  2  to  10  percent,  and  both  symmetrical  and  cam- 
bered profiles.   These  results  are  reported  in  references  2,  3,    and  k   and 
have  been  studied  by  McDevitt  (refs.  3  and  5)  who  showed  in  a  convincing 
manner  that  the  experimental  data  can  be  correlated  successfully  by  means 
of  the  transonic  similarity  rules.   These  same  results  will  be  used  herein 
to  evaluate  one  phase  of  the  transonic  area  rule.   The  relationship  between 
the  two  rules  is  naturally  of  interest  and  will  also  be  explored. 

Since  the  transonic  area  rule  is  considered  to  apply  equally  to  all 
low-aspect-ratio  wing-body  combinations,  detailed  examination  of  such  a 
limited  class  of  aerodynamic  shapes  as  a  family  of  rectangular  wings  is 
not  without  value  inasmuch  as  limitations  revealed  in  special  applications 
must  appear  as  a  limitation  in  the  general  case.  The  restricted  range  of 
the  investigation  is  compensated  somewhat  by  the  fact  that  the  geometric 
simplicity  increases  the  chances  of  understanding  the  underlying  causes. 
Although  the  results  can  only  be  said  to  apply  with  surety  to  the  specific 
cases  investigated,  the  method  of  approach  is  not  restricted  and  may  be 
applied  similarly  to  other  cases  as  more  data  become  available  and  as 
understanding  increases.   In  this  way,  the  present  discussion  may  be  con- 
sidered more  as  suggestive  than  definitive. 


PRINCIPAL  SYMBOLS 


A  aspect  ratio 

A  [(7+1)  M02t]1/3A 


b  wing  span 
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dimensionless  function  describing  thickness  distribution 


g2         dimensionless  function  describing  camber  distribution 
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h  maximum  camber  of  wing 
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M0         Mach  number 

Mref        reference  Mach  number  less  than  the  critical 

q  dynamic  pressure,  -£■  UQ2 
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R  Reynolds  number 

^ref  reference  Reynolds  number  corresponding  to  Mref 

S  area 

Sc  cross-section  area 

Sm  maximum  cross-section  area 

Sp  plan-form  area 

s  dimensionless  area  distribution  function 

t  maximum  thickness 

UQ  free-stream  velocity 

x,y  Cartesian  coordinates  in  plane  of  wing  where  x  extends  in 
the  direction  of  the  free-stream  velocity 

Y  dimensionless  function  describing  plan  form 

ordinates  of  wing  profiles 

Zy  ordinates  of  upper  surface  of  wing  profiles 

Zj  ordinates  of  lower  surface  of  wing  profiles 

a  angle  of  attack 


a 


*o 


a 
t/c 

ratio  of  specific  heats,  for  air  7  =  1.1+ 


Mq2 


[(7  +  1)  M^tF'3 
P0         density 


t 
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DIMENSIONAL  CONSIDERATIONS 


The  forces  on  a  body  moving  through  an  infinite  fluid  are  dependent 
on  a  number  of  parameters  equal  to  that  necessary  to  describe  the  problem. 
Thus,  it  is  well  known  from  dimensional  considerations  that  the  drag  D 
of  a  body  can  be  expressed  as  the  product  of  a  pressure,  say  the  dynamic 

Pn   o 

pressure   q  =  -~   U0  ,  a  characteristic  area  S,  and  some  function  of  Mach 

number  M0,  Reynolds  number  R,  the  geometry  of  the  body,  and  gas  proper- 
ties such  as  the  ratio  of  the  specific  heats  7,  etc.   This  can  be  written 
symbolically  as  follows: 

D  =  q_Sf(M0,  R,  geometry,  gas  properties)  (l) 

(Note  that  S  does  not  necessarily  refer  to  the  plan-form  area,  but, 
rather,  to  any  combination  of  geometric  lengths  having  the  dimensions  of 
area.)   The  geometry  of  a  body  can  be  described  by  a  number  of  dimension- 
less  parameters,  the  number  of  which  depends  on  the  complexity  of  the 
shape.   If  these  considerations  are  applied  to  the  present  family  of  rec- 
tangular wings  having  NACA  63AXXX  profiles,  the  description  is  particularly 
simple  since  the  wing  plan  form  is  determined  by  specifying  the  aspect 
ratio  A,  the  profile  by  the  thickness  ratio  t  =  t/c  and  the  camber  ratio 
H  =  h/c,  and  the  inclination  of  the  wing  by  the  angle  of  attack  a.   In 
this  way,  equation  (l)  can  be  rewritten  in  the  following  more  explicit 
form: 

-g  =  f(M0,  R,  A,  t,  H,  a,  gas  properties)  (2) 

If,  as  in  nearly  all  problems  of  aerodynamic  practice,  all  measurements 
are  made  in  the  same  fluid,  air,  the  gas  properties  can  be  represented  by 
merely  a  set  of  constants  and  therefore  disappear  as  parameters  in  equa- 
tion (2)  leaving  only 


2r  =  f(M0,  R,  A,  r,  H,  a)  (3) 


qS 

If  attention  is  confined  to  the  drag  at  zero  lift  D0,  as  is  the  case  with 
the  transonic  area  rule,  further  simplification  results  because  a  is 
then  a  dependent  parameter  and  may  be  eliminated  even  though  it  is  not  a 
constant  for  variously  cambered  wings.   Then 

2g  =  f(M0,  R,  A,  x,  H)  (1+) 

qS 

The  proof  of  the  foregoing  step  is  as  follows:   The  relationship  for  lift 
analogous  to  equation  (3)  is 

CL  =jg  =   f(Mo,  R,  A,  t,  H,  a)  (5) 
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where  f  is,  of  course,  a  different  function  of  the  indicated  variables. 
This  can  be  solved  for  a 

a  =  f(Mo,  R,  A,  t,  H,  CL)  (6) 

and  substituted  in  its  place  in  equation  (3)  to  produce 

^  =  f(MQ,  R,  A,  T,  H,  CL)  (7) 

In  this  form,  it  is  clear  that  the  restriction  to  zero  lift  fixes  the 
value  of  the  last  parameter,  thereby  leaving  D0/qS  dependent  only  on  the 
remaining  parameters  as  in  equation  (k) . 

Although  dimensional  considerations  such  as  the  foregoing  display  the 
parameters  or  factors  that  influence  the  drag,  they  do  not  provide  any 
information  on  the  nature  of  the  functional  relations  involved.   Thus,  it 
may  be  that  some  of  the  parameters  are  more  important  than  others,  or  that 
a  parameter  is  of  importance  over  a  certain  range  of  values  but  of  negli- 
gible importance  in  another  range,  or  that  a  parameter  is  important  or  not 
depending  on  the  value  of  another  parameter,  etc.   Examples  are  that 
D0/qS  depends  on  Mq  at  transonic  and  supersonic  speeds  but  not  at  sub- 
sonic speeds,  D0/qS  depends  on  A  for  small  A  but  not  for  large  A 
(if  S  refers  to  the  plan-form  area),  etc.   Since  there  exists  a  great 
body  of  literature  pertinent  to  the  nature  of  these  dependencies  which  is 
probably  quite  familiar  to  the  readers,  no  further  amplification  appears 
necessary  here. 

TRANSONIC  AREA  RULE 
Theoretical  Considerations 


The  customary  verbal  statement  of  the  transonic  area  rule  is  given 
in  the  INTRODUCTION.   It  should  be  noted  that  the  term  "drag  rise"  does 
not  refer  to  the  increase  in  drag  force  D0,  but  to  the  increase  in  the 
ratio  D0/q  between  a  subcritical  reference  Mach  number  Mref  and  a 
transonic  Mach  number  Mq.   If  Sc(x)  represents  the  streamwise  or  axial 
distribution  of  cross-section  area,  the  transonic  area  rule  states  that 
the  variations  A(D0/q)  with  Mach  number  MQ  are  the  same  for  all  low- 
aspect-ratio  wing-body  combinations  having  the  same  Sc(x).   Note  that 
the  rule  in  this  form  relates  a  very  wide  class  of  wings,  bodies,  etc., 
but  that  all  related  configurations  must  have  the  same  longitudinal 
length  c.   This  restriction  can  be  removed  easily  by  recasting  the  state- 
ment in  dimensionless  form  using  c  as  the  characteristic  length.   Thus, 
we  can  say  that  the  variation  of  A(D/qc2)  with  Mach  number  Mq  is  the 
same  for  all  low-aspect-ratio  wing-body  combinations  having  the  same 
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Sc(x/c)/c2.   From  a  slightly  different  point  of  view,  the  area  rule  states 
that,  near  the  speed  of  sound,  A(D/qc2)  is  a  function  of  Mq  and 
Sc(x/c)/c2,  or 


r'D0\    ^DQ\     (V0\  [    scU/c) 


(8) 


The  axial  distribution  of  cross-section  area  Sc(x/c)  can  be  written 
in  the  form 

Sc(x/c)  =  Sms(x/c)  (9) 

where  S^  is  the  maximum  cross-section  area  and  s(x/c)  is  an  area- 
distribution  function.  Application  of  the  transonic  area  rule  to  the  drag 
results  of  the  present  family  of  rectangular  wings  is  facilitated  partic- 
ularly by  the  fact  that  s(x/c)  is  the  same  for  all  wings.  As  a  result, 
the  area  distribution  Sc(x/c)  of  any  wing  in  this  family  is  specified  by 
stating  the  value  of  the  maximum  cross-section  area  Sm,  so  that 

A(D0/qc2)  =  f(M0,  Sm/c2)  (lO) 

Since  the  chord  c  and  thickness  t  of  each  wing  are  constant  across  the 
span  b,  the  maximum  cross-section  area  is  equal  to  the  product  of  the 
span  and  thickness 

Sm  =  bt  (11) 

and 

Sm   b  t  ,   , 

q§  =  __  =  AT  (12) 

c    c  c 

whence  equation  (10)  can  be  rewritten  for  this  family  of  airfoils  as 

The  remarkable  simplicity  of  these  statements  is  emphasized  by  com- 
parison with  the  functional  relation  revealed  by  dimensional  consideration 
alone.   Thus,  if  the  unspecified  area  S  is  replaced  by  c2,  equation  (U) 
yields 


A 


(jj*)     =  A(j^>)  =  f(Mo,  R,  A,  t,  H)  -  f(Mref,  Rref,  A,  t,  H)  (lk) 


where  Rref  refers  to  the  value  of  Reynolds  number  associated  with  the 
subcritical  reference  Mach  number  Mref,  and  the  symbol  f  refers  in  each 
case  to  the  appropriate  function  of  the  indicated  variables.   In  the  cus- 
tomary discussion  of  drag-rise  data,  Mref  and  Rref  are  constants  and  no 
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longer  appear  as  parameters  in  equation  (lU).   Further  simplification 
occurs  in  most  cases  because  the  wind-tunnel  or  flight  test  technique 
determines  a  specific  relation  between  the  Mach  and  Reynolds  numbers.  As 
a  result,  either  Mo  or  R  can  be  removed  as  parameters  since  the  value 
of  either  is  determined  by  that  of  the  other.   Since  the  present  problem 
is  more  closely  connected  with  effects  of  compressibility  than  of  vis- 
cosity, it  is  appropriate  to  retain  M0  as  the  significant  parameter.   In 
this  way,  equation  (lh)    reduces  to 


A(qcV   =  f(M°'  A'  T'  H) 


(15) 


Comparison  of  equations  (13)  and  (15)  highlights  the  fact  that  the 
area  rule  affirms,  in  dimensionless  terms,  that  the  drag-rise  parameter 
A(D0/qc2)  for  the  present  family  of  wings  depends  on  Mach  number  and  the 
product  of  aspect  ratio  and  thickness  ratio  At  (or  the  maximum  cross- 
section  area  parameter  Sm/c2)  but  is  independent  of  camber  ratio  H,  and 
aspect  ratio  A  or  thickness  ratio  t  taken  separately. 


Comparison  with  Experiment 


Application  to  wings  having  identical  area  distributions. -  The 
applicability  of  the  transonic  area  rule  to  the  present  family  of  rectangu- 
lar wings  can  be  examined  in  several  ways.   Perhaps  the  most  obvious  way 
is  to  actually  compare  the  variation  of  A(DD/qc2)  with  Mq  for  two  or  more 
wings  having  identical  area  distributions.   For  the  present  family  of 
wings,  this  means  comparing  wings  having  the  same  S™/c2  or  At.   The  tran- 
sonic area  rule  predicts  that  the  variation  of  A(D0/qc2)  with  M0  should 
be  the  same  for  all  such  wings. 

An  example  of  such  a  comparison  is  shown  in  sketches  (a)  and  (b) . 
The  experimental  data  are  from  references  2  and  3-   Both  wings  have 
./Or 


.6  .7  .8 

Sketch  (b) 
CONFIDENTIAL 


NACA  RM  A5UF28 


CONFIDENTIAL 


symmetrical  sections  (H  =  0)  and  At  =  0.l6  t>ut  one  has  an  aspect  ratio  of 
2  and  thickness  ratio  of  0.08;  whereas,  the  other  has  an  aspect  ratio  of 
k   and  thickness  ratio  of  0.0U.  The  first  sketch  shows  the  total  drag  and 
the  second  the  drag  rise  determined  by  subtracting  the  value  of  D0/qc2 
at  Mq  =  0.6.  Although  the  two  curves  in  sketch  (b)  are  not  identical  as 
predicted  by  the  area  rule,  they  are  closely  related.   Innumerable  reasons 
could  be  advanced  to  explain  the  differences  between  the  two  curves;  per- 
haps there  are  viscous  effects  which  may  significantly  affect  the  drag 
rise,  perhaps  the  measurements  are  not  sufficiently  accurate  or  the  flow 
field  sufficiently  uniform,  or  perhaps  the  aspect  ratio  or  thickness  ratio 
is  too  large,  etc.   In  any  case,  this  particular  comparison  would  probably 
be  scored  in  favor  of  the  transonic  area  rule . 


Another  comparison,  this  time  among  three  wings  of  aspect  ratio  2, 
thickness  ratio  0.06,  but  different  amounts  of  camber  is  shown  in  sketches 


(c)  and  (d) 
ficient  Cz. 

designation. 


The  amount  of  camber  is  specified  by  the  ideal  lift  coef- 
in  accordance  with  the  NACA  scheme  for  airfoil  section 
Again,  the  basic  data 


are  presented 
the  drag  rise 


in  the  first  sketch  and 
in  the  second.   In  this 
case,  however,  the  variation  of 
A(D0/qc2)  with  Mq  is  definitely  not 
the  same  for  the  three  wings  and  the 
use  of  the  area  rule  could  lead  to 
serious  error.  At  a  Mach  number  of 
unity,  where  the  area  rule  is  sup- 
posed to  be  most  accurate,  the  drag 
rise  of  the  wing  with  greatest  camber 
is  nearly  twice  that  of  the  uncambered 
wing. 

Wings  having  similar  area  dis- 
tribution.- Although  a  certain  num- 
ber of  additional  comparisons  of  the 
type  described  in  the  preceding  sec- 
tion can  be  made  using  the  data  of 
references  2  through  5>  the  number  is 
definitely  limited  because  the  test 
program  was  not  designed  to  preserve 
a  single  value  for  At  for  all  wings. 
It  is  furthermore  not  practical  to 
carry  out  extensive  programs  of  such 
a  type  because  it  necessitates  the 
testing  of  very  thin  wings  of  high 
aspect  ratio  and  thick  wings  of  low 
aspect  ratio.  As  mentioned  in  the 
derivation  of  equation  (l3),  however, 
all  members  of  the  present  family  of 
wings  have  similar  distributions  of 
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Sketch  (d) 
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cross-section  area  (a  single  area-distribution  function  s(x/c)),  and  the 
transonic  area  rule  can  be  extended  to  include  such  cases  by  introducing 
Sm/c2  or  At  as  a  second  parameter.   In  exchange  for  being  able  to  cor- 
relate the  drag  rise  of  bodies  having  not  only  identical,  but  also  similar 
area  distributions,  we  incur  the  complications  of  a  dependence  on  two 
parameters  rather  than  only  one.   Thus,  whereas  the  curves  representing 
the  variation  of  A(Do/qc2)  with  MQ  for  all  wings  having  identical 
Sjj/c2  or  At  should  coincide  to  form  a  single  line,  those  for  a  family  of 
wings  having  similar  area  distributions  should  form  a  family  of  lines. 
Simplicity  can  be  regained,  however,  by  restricting  attention  to  a  single 
Mach  number  and  ascertaining  the  variation  of  A(D0/qc2)  with  Sm/c2  or  At. 
Since  the  customary  statement  of  the  area  rule  restricts  attention  to 
near-sonic  speeds,  the  most  appropriate  single  Mach  number  to  select  for 
such  a  comparison  is  unity. 

Sketch  (e)  shows  the  variation  of  A(Do/qc2)  with  Sm/c2  for  all  the 
uncambered  wings  of  references  2  and  3«   It  can  be  seen  that  all  these 
results  fall  near  to  a  single  curved  line  for  wings  of  all  aspect  ratios 
up  to  3  hut  that  those  for  wings  of  aspect  ratios  k   and  6  depart  from  this 
line  in  a  systematic  manner. 


Sm/c 


sv  .25 


The  same  results  are  replotted  in  sketch  (f)  versus  the  square  of 
rather  than  the  first  power.   It  can  be  seen  that  the  curved  line 
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of  sketch  (e)  for  wings  having  aspect  ratios  less  than  3  is  now  a  straight 
line,  indicating  that  the  sonic  drag  rise  for  these  wings  is  proportional 
to  the  square  of  the  maximum  cross-section  area.   Such  a  dependence  is 
consistent  with  the  formulas  of  linearized  compressible  flow  theory  for 
the  wave  drag  of  nonlifting  slender  wings  and  bodies  at  supersonic  speeds. 

These  plots  can  be  interpreted  as  showing  that  the  sonic  drag  rise 
of  the  uncambered  members  of  the  present  family  of  wings  depends  on  the 
cross-section  area  in  accordance  with  the  area  rule,  provided  the  aspect 
ratio  is  about  3  or  less.   On  the  other  hand,  the  results  for  wings  of 
higher  aspect  ratio  can  only  be  interpreted  as  indicating  that  some  parame- 
ter other  than  the  cross-section  area  must  be  involved.   Inasmuch  as  the 
geometry  of  wings  of  the  present  family  is  described  completely  by  the  two 
parameters,  aspect  ratio  and  thickness  ratio,  it  is  clear  that  these 
parameters  must  assume  importance  in  some  form  other  than  their  simple 
product  At  for  wings  of  larger  aspect  ratio.   One  could  seek  the  new 
relation  empirically,  but  the  transonic  similarity  rules  provide  a  theo- 
retical basis  for  proceeding.   Some  properties  of  transonic  similarity 
rules  are  reviewed  in  the  following  section,  although  the  reader  is  refer- 
red to  the  original  references  for  further  details. 


TRANSONIC  SIMILARITY  RULES 
Statement  of  Rule 


Transonic  similarity  rules  are  derived  from  the  nonlinear  equation 
of  inviscid  flow  theory  and  are  known  for  thin  wings  (e.g.,  ref .  6)  and 
slender  bodies  of  revolution  (ref.  7)>  but  not  for  wing-body  combinations 
with  pointed  noses.   In  contrast  to  the  transonic  area  rule  which  relates 
the  zero-lift  drag  rise  of  families  of  bodies  having  identical  or  similar 
axial  distributions  of  cross-section  area,  the  transonic  similarity  rules 
relate  the  aerodynamic  properties  of  much  more  highly  restricted  families 
of  bodies.  Even  for  wings  alone,  the  restrictions  imposed  on  the  members 
of  a  single  family  are  much  more  severe  than  for  the  area  rule,  since  all 
members  of  a  single  family  must  have  affinely  related  plan  forms,  affinely 
related  thickness  distributions,  and  affinely  related  camber  distributions. 
To  be  more  explicit:   if  the  plan  form  is  given  by 


^2  =  Y(x/c)  <16) 
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as  indicated  in  sketch  (g),  it  is  required  that 
Y(x/c)  be  a  single  function  for  all  wings  of  a 
given  family.   Furthermore,  if  the  ordinates  of 
the  wing  surface  are  given  by 


2  =  *  Li 

c   c  ~2 


*  1 


h 
+  t 


1     I 


a  x 
t/c  c 


(IT) 


Sketch  (g) 


where  the  plus  sign  is  associated  with  the  ordinates  of  the  upper  surface 
Zu  and  the  minus  sign  with  those  of  the  lower  surface  Zz,  it  is  required 
that  g1(x/c,  y/b)  and  g2(x/c,  y/b)  be  single  functions  for  all  wings  of 
a  given  family.   The  first  of  these  restrictions  requires  that  related 
plan  forms  be  obtainable  one  from  another  by  a  differential  lengthening  of 
lateral  and  longitudinal  dimensions.   Thus  rectangular  plan  forms  consti- 
tute one  family,  triangular  plan  forms  with  straight  trailing  edges 
another,  etc.   Examples  of  related  plan  forms  are  illustrated  in 
sketch  (h) .   The  relationship  expressed  in  equation  (l7)  requires  that  the 

thickness  distribution  and  camber 


□ 


-> 


0 

Sketch  (h) 


/\    variations  must  be  the  same  for  all 
wings  of  a  particular  family.   The 
magnitudes  of  the  maximum  thickness 

A    ratio  t/c  and  camber  thickness 
ratio  h/c,  as  well  as  the  angle  of 


attack 


a, 


may  be  different  for  vari- 


ous members  of  a  single  family. 


From  the  present  point  of  view,  one  of  the  most  significant  proper- 
ties of  the  family  of  wings  described  above  is  that  all  members  have  the 
same  dimensionless  area-distribution  function  s(x/c).   This  can  be  shown 
as  follows: 


+Y(x/c)  +Y 

S(x/c)  =  J  (Zu  -  Zz)dy  =  I 

J-Y(x/c)  J-Y 


tfc  (§,*)* 


=  bt 


+Y/b 

r 

-Y/b 


f'f 


=  btg(x/c)  =  SmsU/c)    (18) 
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where  g  represents  the  integral  of  the  preceding  expression  and  is  pro- 
portional to  s(x/c),  since  Sm  is  proportional  to  bt  for  a  family  of 
wings.  Because  of  this  fact,  a  family  of  wings  suitable  for  correlation 
by  the  transonic  similarity  rules  is  also  part  of  a  family  of  bodies  suit- 
able for  correlation  by  the  area  rule. 

For  wings  of  such  a  family,  the  similarity  rules  of  inviscid,  small- 
disturbance,  transonic  flow  theory  provide  that  the  wave  drag  LV  be 
given  by 

T5/3 
=    ^P     [Mo%+l)]l/3    f  (t0>     K>    £'     S)  U9) 

where 

T  t/c 

Mo2 


&0        [(7+l)M02T]2/3 

A        [(7+l)MoS-]1/&A 

h       h/t 

a       a/r 

and  S„  is  the  plan-form  area  of  the  wing.   Equation  (19)  is  not  in  the 
proper  form  for  comparison  with  the  area  rule  since  the  latter  is  con- 
cerned with  conditions  at  zero  lift  and,  hence,  imposes  an  indirect 
requirement  on  a.  We  can  proceed  toward  the  desired  form,  however,  by 
introducing  the  transonic  similarity  rule  for  the  lift  coefficient 
CL  =  L/qSp 

°L  =  [M^l)]^  f(*o'  l>   *>   °>  (20) 


and  defining  the  reduced  lift  coefficient 

[M0g(y+l)]1/3 

T2/3 


CL  -  L"°  "+i"    CL  =  f  (*  ,  A,  h,  a)  (21) 


Now  since  C^  is  a  function  of  the  same  four  variables  that  appear  in 
equation  (l9)>  cc  can  be  replaced  with  Cl  i-n  the  latter  equation,  whence 

°w  =  *Sp  [^g^Jl/a  f  (t0,  A,  h,  CL)  (22) 
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The  condition  of  zero  lift  eliminates  the  last  parameter,  leaving  only 


D 


.5/3 


Ov 


=  qSp  i^Li)]!'*  f{to>  l> 5) 


qS 


-5/3      r   Mq2 


P    [Mo2(7+Dl1/3  [[(7+-)Mo2t]2/ 

[(7+l)Mo2T]1/3A,    h/t! 


(23) 


Application  to  Rectangular  Wings  and  Comparison 
with  Transonic  Area  Rule 


The  restrictions  introduced  in  the  derivation  of  the  transonic  simi- 
larity rule  for  zero-lift  wave  drag  are  such  that  they  permit  the  direct 
application  of  equation  (23)  to  the  present  family  of  rectangular  wings 
having  NACA  63AXXX  profiles.   It  is  natural  to  compare  this  functional 
relation  with  the  corresponding  relations  of  equations  (lO)  and  (l3)  given 
by  the  transonic  area  rule. 


AG?) =  f(M°' Sm/c2)  =  f(M°'  at) 


(2k) 


At  first  glance,  the  two  sets  of  relationships  appear  to  bear  only  slight 
resemblance.   It  can  be  seen  upon  closer  examination,  however,  that  some 
of  the  apparent  differences  are  superficial  and  of  little  or  no  signifi- 
cance.  For  instance,  equation  (23)  is  concerned  with  wave  drag  Dq^. 
whereas  equation  (2k)    is  concerned  with  drag  rise.   It  is  evident,  how- 
ever, that  the  two  rules  are  actually  concerned  with  the  same  quantity, 
since  the  drag  rise  can  be  considered  to  be  an  approximation  for  the  wave 
drag  under  the  assumption  that  the  friction  drag  coefficient  Is  independent 
of  Mach  number. 

Another  point  of  apparent  lack  of  resemblance  is  that  equation  (23) 
does  not  show  an  explicit  relationship  between  wave  drag  and  maximum 
cross-section  area,  since  the  latter  is  not  a  function  of  SpT5^3  nor  of 
|0,  A,  or  h  taken  separately.   This  does  not  mean  that  the  transonic 
similarity  rule  is  incompatible  with  the  area  rule  because  there  are  sev- 
eral possibilities  for  making  such  a  dependence  visible.   Two  permissible 
procedures  are  to  multiply  the  right  side  of  equation  (23)  by  either  iQ 
or  A.   The  first  procedure  is  of  no  help  for  the  discussion  of  conditions 
at  M0  =  1,  since  the  quantity  M02  -  1  appears  in  two  places  and  an  inde- 
terminate form  ensues.  The  second  is  perfectly  acceptable,  however,  and 
produces  the  following  relationships: 
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Dow  =  qSpAr2  f(|Q,  A,  h) 


qC\^y)   f(^0'  K>    5)  (25) 


Since  the  maximum  cross-section  area  Sjq  of  the  present  family  of  rec- 
tangular wings  is  equal  to  the  product  bt,  equation  (25)  can  be  rewritten 
as  follows: 

Dov  =  ^c2(^i)2  fU0,  A,  h) 


c2/      1[(7+i)Mo2t]2/3' 

[(7+1)Mo2t]1/3  A,   h/tl  (26) 


Equation  (2*0  can  likewise  be  rewritten  by  multiplying  f  by  the  square 
of  Sm/c2, 


A 


&)  "  (¥J   f(M°'  "^  '  (?J  »(*.  *T)        (27) 


This  appears  to  be  the  closest  that  the  two  rules  can  be  brought 
together  without  introducing  additional  simplifications.  Both  rules  are 
now  concerned  with  essentially  the  same  quantity,  D^/qc2  and  A(D0/qc2)  . 
Each  rule  states  that  this  quantity  is  proportional  to  the  square  of 
Sm/c2  times  some  unknown  function  of  certain  specified  parameters.   The 
two  rules  disagree  completely,  however,  as  to  the  nature  of  the  parameters. 
The  transonic  similarity  rule  specifies  three  parameters  iQ,   A,  and  h, 
whereas  the  area  rule  specifies  only  two,  Mq  and  At.   Since  neither  set 
can  be  transformed  into  the  other,  it  is  apparent  that  the  only  way  in 
which  both  rules  can  be  universally  correct  is  for  the  function  f  to  be 
actually  a  constant  and  not  to  depend  on  the  value  of  any  of  the  five 
parameters.   It  is  obvious  that  such  is  not  the  case,  however,  since  it 
requires,  for  instance,  that  the  wave  drag  be  independent  of  Mach  number. 

There  is  a  way  out  of  this  apparent  impasse,  however,  if  the  range 
of  validity  of  one  or  both  of  the  rules  is  restricted  sufficiently  that 
f  is  independent  of  the  remaining  variables.  Comparison  of  the  two  rules 
shows  that  the  drag  depends  in  both  cases  on  the  Mach  number  of  the  stream 
and  the  geometry  of  the  wing  and,  in  the  case  of  the  transonic  similarity 
rule,  on  the  ratio  of  the  specific  heats  7.   Since  nearly  all  problems  of 
aerodynamic  interest  are  concerned  solely  with  air,  however,  y      is  a  con- 
stant and  need  not  be  retained  as  a  parameter.   The  resulting  simplified 
equation  can  be  written  in  full  in  either  of  the  following  forms: 
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r  »/r  2 


£ss .  ra»Y  f  r,»i^.  (vm^  *,  vt 


2    V   2 


2T\2/3; 


=  (AT)2  f 


(Mq2T) 

Mo2  -  1    ,„  ,  ,1/3      / 
(„  2T)2/3^  (Mo  T)     A,  h/t 


(28) 


With  regard  to  further  simplifications,  it  should  be  noted  that  the 
derivation  of  the  transonic  similarity  rule  requires  that  the  thickness 
and  camber  be  small  with  respect  to  the  chord,  but  does  not  restrict  the 
Mach  number  or  aspect  ratio.   The  transonic  area  rule,  on  the  other  hand, 
restricts  itself  to  Mach  numbers  near  the  speed  of  sound  and  to  wings 
having  low  aspect  ratios. 

The  functional  relation  given  in  equation  (27)  representing  the  area 
rule  simplifies  if  attention  is  fixed  on  any  given  Mach  number,  since  Mq 
appears  as  an  isolated  parameter,  thus 

Do  \    /Sm 


On  the  other  hand,  the  only  Mach  number  at  which  the  corresponding  rela- 
tion given  by  the  transonic  similarity  rule  simplifies,  irrespective  of 
the  thickness  ratio,  is  Mq  =  1.   The  resulting  expression  is 

22g  =  (J|>  f(AT1/3,  h/t)  (30) 

This  expression  simplifies  further  for  wings  having  symmetrical  sections, 
since  h/t  is  zero.   For  such  cases,  the  transonic  similarity  rule  for 
wave  drag  at  Mq  =  1  reduces  to 

ASmY  f(ATi/3)  (31) 


qc2    WJ 

The  foregoing  analysis  has  developed  certain  equations  relating  to 
the  drag  rise  or  wave  drag  of  a  family  of  wings  having  relatively  simple 
geometry.   Despite  the  great  restrictions  imposed  by  the  selection  of  such 
a  family  of  bodies,  the  analysis  has  disclosed  a  number  of  significant 
points  which  complement  those  discussed  in  connection  with  the  experimental 
data  presented  in  sketches  (a)  through  (f).   First  of  all,  it  has  been 
shown  that  the  only  way  in  which  the  wave  drag  at  Mq  =  1  can  depend  on 
the  axial  distribution  of  cross-section  area  (defined  in  the  present  fam- 
ily of  wings  by  the  value  of  the  maximum  cross-section  area  Sm)  and  still 
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be  compatible  with  the  transonic  similarity  rule  for  drag  is  for  Dov/qc2 
to  be  proportional  to  the  square  of  Sm/c2.   This  is  precisely  the  depend- 
ence disclosed  by  the  experimental  data  for  low-aspect-ratio  wings  having 
symmetrical  sections  and  shown  in  sketch  (f).   The  transonic  similarity 
rule  states  that  the  wave  drag  may  depend  on  the  camber,  whereas  the  area 
rule  states  that  it  does  not.   The  experimental  data  in  sketch  (d)  show 
that  camber  has  a  significant  effect  on  drag. 

It  is  evident,  both  from  a  priori  considerations  and  from  the  experi- 
mental results  shown  in  sketch  (f),  that  some  change  must  occur  in  the 
relation  between  wave  drag  and  maximum  cross-section  area  as  the  aspect 
ratio  becomes  very  large.  The  only  possibility  permitted  by  the  transonic 
area  rule  is  that  A(D0/qc2)  varies  with  Sm/c2  in  some  other  manner  than 
as  the  square.   It  can  be  seen  from  sketch  (e),  however,  that  the  data  for 
wings  of  the  present  family  cannot  be  correlated  on  this  basis  if  the 
aspect  ratio  is  greater  than  about  3-   The  transonic  similarity  rule  for 
the  zero-lift  wave  drag  of  uncambered  wings  at  Mq  =  1  provides  a  differ- 
ent dependence  by  stating  that  Dc^/qc2  is  equal  to  the  square  of  Sm/c2 
times  some  function  of  At1/3.   This  statement  is  compatible  with  the 
area  rule  if  Dov/qc2  is  independent  of  At1/3  for  small  values  of  the 
latter.   If  wave-drag  data  for  wings  of  larger  aspect  ratio  can  be  cor- 
related successfully  by  considering  the  value  of  this  quantity,  a  theoret- 
ical basis  for  a  limit  to  the  range  of  applicability  of  the  transonic  area 
rule  has  been  found. 


Comparison  with  Experiment 


The  foregoing  discussion  has  served  to  focus  attention  on  the  fact 
that  the  parameter  At  /3,  already  familiar  from  prior  papers  on  transonic 
flow  (e.g.,  ref.  3)>  may  be  of  importance  in  defining  the  limit  of  appli- 
cability of  the  transonic  area  rule  as  applied  to  a  family  of  af finely 
related  wings.   The  functional  relation  of  equation  (31)  suggests  that  if 
D^c^qSm2  is  plotted  as  a  function  of  At1/3,  all  the  data  for  the 
uncambered  wings  should  fall  on  a  single  curve.  Moreover,  the  values  of 
D(v. c2/qSm2  should  be  independent  of  At1 /3   over  whatever  range  the  area 
rule  applies.  Although  this  method  of  plotting  is  conceptually  simple, 
it  imposes  severe  requirements  on  the  accuracy  of  the  experimental  deter- 
mination of  wave  drag  because  any  errors  are  magnified  as  a  result  of 
dividing  by  the  square  of  the  maximum  cross -section  area.   For  results 
such  as  the  present  in  which  the  wave  drag  is  not  actually  measured,  but 
is  inferred  from  measurements  of  total  drag  by  subtracting  an  estimated 
friction  drag,  greatest  difficulties  are  experienced  if  the  wings  are  thin 
and  of  low  aspect  ratio  so  that  the  wave  drag  is  only  a  small  fraction  of 
the  total  drag.   Consequently,  data  points  for  which  the  wave  drag  is  less 
than  half  the  total  drag  are  omitted  in  the  plot  of  A(D0c2/qSm2)  versus 
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AT1 'a  shown  in  8ketch  (i).   Even  with 
this  precaution,  the  data  evidences 
considerable  scatter  for  thin,  low- 
aspect-ratio  wings.   (The  symbols  refer 
to  the  same  wings  as  in  sketch  (e).) 

The  principal  points  of  interest 
in  this  plot  are  threefold.   First, 
except  for  the  smallest  values  of 
At  3  where  the  scatter  is  too  large 
to  provide  any  positive  conclusions, 
the  points  determine  essentially  a 
single  line,  indicating  that  the  sonic 
wave -drag  characteristics  of  the  uncam- 
bered  wings  of  the  present  family  can 
be  correlated  successfully  by  the  tran- 
sonic similarity  rule  (this  has  been 
shown  previously  by  McDevitt,  ref.  3);  second,  A(D0c2/qSm2)  is,  at  best, 
independent  of  At1/3  for  values  of  the  latter  up  to  about  unity,  indi- 
cating that  the  drag  in  this  range  varied  in  accordance  with  the  area  rule 
as  well;  third,  A(D0c2/qSm2)  varies  appreciably  with  At1/3  at  larger 
values  of  the  latter,  indicating  that  the  range  of  applicability  of  the 
area  rule  was  exceeded. 
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Sketch  (i) 


Before  leaving  the  discussion  of  the  correlation  of  drag  data  at 
Mq  =  1  by  use  of  the  transonic  similarity  rule,  it  is  worthwhile  to  call 
attention  to  another  method  of  plotting  which  illustrates  the  same  factors 
although  in  a  slightly  different  manner.   This  procedure,  which  is 
described  more  completely  in  references  3>  5>  and  6,  is  based  on  equa- 
tion (23)  rather  than  equation  (26)  and  consists,  for  the  zero-lift  drag 
of  a  family  of  uncambered  wings  in  an  air  stream  with  Mq  =  1,  of  plotting 
the  variation  of  Cj)~ /t5/3  with  At1/3  where  Cj)~  rePresen"ts  ^W/qSp. 

Sketch  (j)  shows  the  data  of  sketch  (i) 
replotted  in  this  manner.   (Once  again, 
the  lack  of  wave-drag  data  requires  the 
substitution  of  drag-rise  information 
and  the  introduction  of  the  drag-rise 
coefficient  ^Dn>    defined  as  equal  to 
A(D0/qSp).) 

The  curve  formed  by  the  data 
points  on  this  type  of  plot  is  perhaps 
somewhat  simpler  than  that  of  sketch 
(i)  since  it  is  asymptotic  to  straight 
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lines  at  both  large  and  small  At 
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the  points  define  a 


Sketch  (j) 


straight  line  passing  through  the 
origin.   Such  a  line  is  in  accordance 
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with  the  area  rule.  At  values  of  At1/3  larger  than  about  unity,  how- 
ever, the  line  determined  by  the  data  points  departs  from  this  initial 
trend  and  turns  toward  the  horizontal.   This  trend  is  contrary  to  the 
area  rule  but  consistent  with  the  fact  that  the  results  for  wings  of  high 
aspect  ratio  must  tend  toward  those  for  wings  of  infinite  aspect  ratio. 
It  cannot  be  emphasized  too  much  that  the  critical  value  of  unity  for 
At1/3  is  determined  solely  on  the  basis  of  data  for  a  very  special  family 
of  rectangular  wings  having  symmetrical  profiles.   Other  families  of  wings 
would  be  represented  by  different  curves  on  such  a  plot. 


It  is  interesting  to  consider  for  a  moment  the  nature  of  this  limit 
and  to  compare  it  with  the  verbal  restriction  of  the  area  rule  to  low- 
aspect-ratio  wings.  Although  the  one  requires  that  AtI^3  t>e  small,  and 
the  other  that  A  be  small,  these  two  statements  are  in  better  agreement, 
insofar  as  engineering  applications  are  concerned,  than  might  appear  at 
first  glance.   This  agreement  results  from  the  fact  that  other  consider- 
tions,  such  as  designing  for  structural  strength  or  for  the  avoidance  of 
excessive  drag,  tend  to  preserve  a  rather  narrow  range  for  the  values  of 
t  likely  to  be  met  in  practice.   The  effectiveness  of  t  as  a  parameter 
is  further  diminished  by  the  fact  that  only  its  cube  root  is  involved. 
Consequently,  the  restriction  to  small  At*'3  represents  a  limitation 
primarily  on  the  aspect  ratio  and  only  secondarily  on  the  thickness  ratio. 


A=I,T  =  04,  Ar*=342     A  =2,  r'./O,  At*=  928 


It  is  also  interesting  to  exam- 
ine the  shapes  of  equivalent  bodies  of 
revolution  having  the  same  axial  dis- 
tribution of  cross-section  area  as 
wings  lying  on  either  side  of  the 
limit.  Accordingly,  sketch  (k)  has 
been  prepared  showing  the  shapes  of 
bodies  of  revolution  equivalent  to 
one  of  the  wings  having  At1/3  much 
less  than  unity,  to  the  two  wings 
having  At1 /3  nearest  unity,  and  to 
the  wing  having  the  largest  At1/3  of 
any  tested.   It  can  be  seen  that  the 
equivalent  bodies  are  blunt  and  stubby 
rather  than  pointed  and  slender.   Thus, 
although  it  has  been  shown  that  the 
drag-rise  characteristics  of  the  var- 
ious members  of  the  present  family  of  uncambered  rectangular  wings  are 
related  to  one  another  in  the  manner  predicted  by  the  transonic  area  rule, 
provided  At1  3  is  less  than  about  unity,  it  would  appear  that  the  drag- 
rise  characteristics  of  the  equivalent  bodies  of  revolution  might  be  con- 
siderably different. 

It  is  shown  in  sketch  (d)  that  camber  has  an  effect  on  the  zero-lift 
drag  rise.   The  transonic  similarity  rule  suggests  that  the  ratio  of  maxi- 
mum camber  to  maximum  thickness  h/t  might  be  an  appropriate  parameter  to 


A=3,r=  04,  Ar*= 1.026      A=6,r-  10,  AT*=2784 

Sketch  (k) 
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use  in  addition  to  At1/3  to  correlate  the  sonic  zero-lift  wave-drag 
characteristics  of  a  family  of  cambered  wings.   Accordingly,  sketch  (z) 

has  been  prepared  showing  the  influence 
of  h/t  for  wings  of  various  At1 '3. 
For  reasons  of  simplicity,  only  three 
values  for  h/t,  namely,  0,  0.222,  and 
O.Mj-U,  are  included  on  this  plot.   The 
important  effects  of  camber  are  readily 
evident  from  this  graph.   It  is  appar- 
ent that  the  area  rule  is  not  applica- 
ble to  wings  having  different  camber- 
thickness  ratios,  even  if  the  values  of 
At  3  are  sufficiently  small  to  permit 
successful  correlation  of  the  wave  drag 
of  uncambered  wings.   On  the  other  hand, 
it  is  permissible  in  the  theory  and 
seems  to  be  indicated  by  the  experimen- 
tal data  that  the  area  rule  is  applica- 
ble to  families  of  cambered  wings  pro- 
viding h/t  is  maintained  constant. 
This  result  is  recognized  in  sketch  (Z) 
by  the  fact  that  A(CD  /t5/3)  for  the 
wings  of  constant  h/t  is  approximately 
proportional  to  At1/3. 
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SUMMARY  OF  RESULTS 


The  range  of  applicability  of  the  transonic  area  rule  has  been  inves- 
tigated by  comparison  with  the  appropriate  similarity  rule  of  transonic 
flow  theory  and  with  available  experimental  data  for  a  large  family  of 
rectangular  wings  having  NACA  63AXXX  profiles.   These  wings  are  of  af finely 
related  geometry  and  are  hence  immediately  amenable  to  analysis  by  the 
transonic  similarity  rules.   On  the  other  hand,  the  axial  distributions  of 
cross-section  area  are  not  identical,  in  most  cases,  but  merely  similar. 
(The  ratio  of  the  local  cross  section  to  the  maximum  cross  section  is  a 
given  function.)   It  is  shown,  however,  how  the  transonic  area  rule  can 
also  be  used  to  correlate  the  sonic  drag-rise  data  for  such  a  family  of 
wings. 


It  is  found  that  the  sonic  zero-lift  drag-rise  data  for  the  present 
family  of  wings  can  be  successfully  correlated  on  the  basis  of  the  area 
rule,  provided  the  wing  profiles  are  symmetrical  and  the  product  of  the 
aspect  ratio  and  the  cube  root  of  the  thickness  ratio  is  less  than  about 
unity.   Within  this  range,  the  sonic  drag  rise  varied  as  the  square  of  the 
maximum  cross-section  area,  all  wings  having  equal  chords.   It  is  demon- 
strated that  this  is  the  only  dependence  of  drag  on  maximum  cross-section 
area  for  a  family  of  wings  like  the  present  that  is  compatible  with  both 
the  area  rule  and  the  transonic  similarity  rule. 

CONFIDENTIAL 


NACA  RM  A5UF28  CONFIDENTIAL  21 


It  was  found  that  the  addition  of  camber  greatly  increased  the  sonic 
drag  rise  and  that  the  application  of  the  transonic  area  rule  to  a  family 
of  wings,  some  of  which  are  cambered  and  others  not,  could  lead  to  serious 
error.   On  the  other  hand,  it  is  indicated  by  the  transonic  similarity 
rules  and  the  experimental  data  that  the  area  rule  is  applicable  to  fami- 
lies of  cambered  wings,  provided  the  camber  distribution,  as  well  as  the 
area  distribution,  are  similar  and  that  the  ratio  of  the  maximum  ordinates 
of  the  camber  and  thickness  distribution  is  maintained  constant. 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  June  28,  1951*- 
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